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1 . INTRODUCTION 


This report is intended as a guide to the users of the inner radia- 
tion zone electron model AE-5. Section 2 includes a description of the 
model, the forms in which it is available, directions on how to use the 
model, and a discussion of its limitations. Computer programs MODEL 
and ORP are described in Sections 3 and 4, respectively. These are 
major programs needed to use the electron models AE-4 and AE-5 and the 
smoothed proton models . 

This document is a companion to one published previously by Teague 
and Vette ("The Inner Zone Electron Model AE-5," NSSDC 72-10, 1972), 
and the reader is referred to that document for a complete description 
of the development of the model. Work is currently in progress to 
improve the high-energy part of model AE-5. In addition, a new proton 
model is being developed. When both of these models are completed, 
they will be compatible with the computer programs described in this 
document . 
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2. DESCRIPTION, USE, AND LIMITATIONS OF MODEL AE-S 


A. Description of the Model 

Model AE-5 describes the inner radiation zone electron environment 
and is based on data from five satellites spanning the period December 
1964 to December 1967. The model provides omnidirectional integral 
flux for energy thresholds in the range 4.0 > Ej/(MeV) > 0.04 and 
for L values in the range 2.8 > L/(R e ) * 1.2 for an epoch of October 
1967. Confidence codes for certain regions of B-L space and certain 
energies are given based on data coverage and the assumptions made in 
the analysis. 

Data from satellites 0G0 1, 0G0 3, 1963 -38C, 0V3-3, and Explorer 
26 were used. The University of Minnesota electron spectrometers were 
carried on board OGOs 1 and 3 and produced data used for model AE-5 
(NSSDC data sets 64-054A-21A and 66-049A-22A supplied to the Data 
Center by Prof. John Winckler and Dr. Karl Pfitzer). These measure- 
ments extended over the period September 1964 to December 1967 . The 
1963-38C satellite was launched in September 1963 and provided data 
through 1967. This spacecraft carried an integral electron spectrometer 
from the Applied Physics Laboratory (Beall, 1969). Data obtained from 
mid-1966 to late 1967 were used in developing model AE-5. Explorer 26 
data from detectors designed by Mcllwain were used for the time inter- 
val January to June 1965. 0V3-3 data from the Aerospace Corporation 

differential nine-channel electron spectrometer, supplied by Vampola 
late in the development of the model, were also incorporated in model 
AE-S. 


The model forms available to a user include a graphical presenta- 
tion and a variety of computer programs. This section describes the 
types of graphs and includes examples of each type. Computer programs 
are discussed in Sections 3 and 4. 
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In previous documentation on trapped particle models, the major 
display has been in the form of omnidirectional integral flux tables. 
Model AE-5 , however, is presented in the form of two-dimensional carpet 
plots J = J(B,L) for given energy thresholds. In addition, carpet plots 
are used for the graphical presentation of the solar cycle variation 
expressed as a ratio. While the omnidirectional flux data have been 
presented tabularly in previous model documentation with greater reso- 
lution than can be obtained from carpet plots, the error associated with 
determining a number from the carpet plots is considered insignificant 
in comparison to the inherent error associated with the model. 

Omnidirectional flux plots are presented in Figures 1 through 7 
for threshold energies Ex = 40, 100, 250, and 500 keV and 1, 2, and 4 
MeV. Fluxes at nongrid B, L, and E points may be obtained simply by 
interpolation as described in Appendix A of Teague and Vette (1972). 

B-L and R-A flux maps are presented in Figures 8 through 13 for thresh- 
old energies Ex = 40 keV, 500 keV, and 1 MeV. In addition, a physical 
impression of the model at these energies may be obtained from the 
three-dimensional plots given in Figures 14 through 16. While the basic 
epoch of model AE-5 is October 1967 corresponding approximately to solar 
maximum, AE-5 contains approximate values of the solar cycle parameter 
for time T 


R t (E t ,L,T) 


_ J(E T ,L,T = October 1967) 
J(E T ,L,T) 


Plots of Rx(L,T) are presented in Figures 17 through 20 for energy 
thresholds Ex = 40, 100, 250, and 500 keV. In these plots, the time 
parameter T has units of months from solar minimum taken as September 
1964. Values of Rx are not presented for higher energies because of 
magnetic storm effects (Teague and Vette, 1972). While the values of 
Rx presented in Figures 17 through 20 have been determined from data 
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over the period 1964 to 1967, they may be used to obtain very approximate 
estimates of the solar cycle effects for epochs later than October 1967. 
Using as a basis the Zurich Sunspot Number, it may be assumed that the 
flux is constant until T = 69 approximately (June 1970) and thereafter 
decreases approximately as described by Figures 17 through 20, reaching 
a minimum at approximately T * 100 (January 1973). It should be ap- 
preciated that extrapolating the model solar cycle dependence in this 
manner is likely to provide very approximate flux estimates only (Sec- 
tion 2B) . 

B. Limitations Of the Model 

It should be remembered that model AE-5 is presented for an epoch 
of October 1967, and temporal variations may result in significant flux 
changes in certain regions of B-L-E space. These temporal variations 
include magnetic storm effects, solar cycle effects, and the decay of 
residual Starfish electrons. These effects are discussed in detail in 
the paper by Teague and Vette (1972). With the exception of magnetic 
storm effects, these temporal variations, however, cause the flux to 
decrease from that given by AE-5 at epoch October 1967, and thus the 
model provides a conservative estimate of the influence of trapped inner 
zone electrons on orbital vehicles. 

To enable the user to assess the reliability of model AE-5, a sys- 
tem of confidence codes is presented. In developing these codes a num- 
ber of criteria were used: number of data sets used, data coverage, 

the degree of data agreement, errors introduced by modeling technique, 
and uncertainties introduced by temporal variations. A scale of 1 to 
10 is used, where 10 corresponds to the highest reliability with an 
expected error of 2 or less and 1 corresponds to the least reliability 
with an expected error in excess of a factor of 10. In general, how- 
ever, efforts have been made to provide pessimistic flux estimates where 
low confidence codes are given so that it is more probable that 
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the flux is lower than the quoted value than higher. Two sets of codes 
are given -- one for the omnidirectional flux at an epoch of October 
1967 (Table 1) and one for the solar cycle parameters (Table 2) . In 
each case, a brief explanation for the confidence code and a section 
reference to the paper by Teague and Vette (1972) are given. 

For example, in Table 1, where the B range is >Bo, the L range 
is 1.2 to 1.4 Rej and the Ej range is >3 MeV, the confidence code is 
1. This code indicates that in these ranges the omnidirectional flux 
indicated by model AE-5 has low reliability, with an expected error 
of more than a factor of 10. The comment column indicates that this 
error results from extrapolation on both B dependence and spectrum 
and from a lack of data. 

Three temporal variations have been noted in the inner radiation 
zone: the decay of Starfish electrons, solar cycle effect, and magnetic 
storm effects. AE-5 attempts to model all three of these, and the 
reader is referred to Teague and Vette (1972) for a complete description 
of the modeling techniques used. 

Model AE-5 contains a small Starfish residual flux in the energy 
range 500 keV < E-p 5 3 MeV and the L range 1.2 < L/ (R e ) < 1.5. Because 
of the lack of Starfish-free data, natural flux levels could not be 
obtained in these intervals, and corresponding low confidence limits 
are quoted in Table 1. Estimates have been made of the times at which 
the Starfish flux component has decayed to the level of the natural 
flux component, and these axe presented in Table 3. For the L and Ej 
region of AE-5 that is influenced by Starfish electrons, it is estimated 
that a maximum reduction of a factor of 5 will result from the decay of 
this component. 
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Some discussion of the solar cycle effect has been given in the 
previous section. The confidence limits presented in Table 1 are ap- 
plicable to the model at epoch October 1967. At other epochs the 
confidence codes are smaller because of the solar cycle effect. If 
Figures 17 through 20 are used to estimate this solar cycle effect, 
however, only the higher confidence codes should be reduced. That is, 
where the model is already associated with a factor of 5 or 6 error, 
no further error is introduced, whereas errors of a factor of 2 will be 
increased to 3 or 4 dependent upon the value of R-p. 

The third temporal variation included in model AE-5 is the effect 
of magnetic storms. This effect is most noticeable at L = 1.9 to 
2.4 R e , and Ej = 0.4 to 2 MeV. Three variables are considered in 
determining the magnetic storm effect: the frequency of the storm, 

the intensity in relation to the undisturbed (quiet day) background, 
and the duration. Assessment of the importance of magnetic storms can 
be performed in practice with consideration of the first two variables 
alone because these exhibit much greater variation with E and L than 
does the third variable. 

The frequency of magnetic storms in the inner belt is too low for 
a statistical approach. However, although the storms are infrequent, 
their relative intensity is high. In these circumstances the flux 
varies considerably from quiet to storm conditions in such a way that 
the changes from one condition to another are unpredictable and cannot 
easily be modeled. 

A crude method has been adopted in model AE-5 for including mag- 
netic storm effects. Average fluxes including magnetically disturbed 
and quiet periods were determined for the period June 1966 to December 
1967. Magnetic storm effects were found to influence this average in 
the region L > 1.8 and 0.4 < Ej/(MeV) < 2.0. The maximum ratio of 
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average to quiet period flux was found to be 40 approximately at 
Ej = 1 MeV at L= 2.4 R e . At L = 1.9 R^, this ratio had reduced 
to 3. It should be appreciated that the average storm flux included 
in AE-S provides inaccurate estimates of the instantaneous fluxes and 
an inaccurate basis for orbit flux integrations because of the low 
frequency and high intensity of magnetic storm effects. These inac- ' 
curacies are incorporated in the confidence limits presented in Table 
1. A model is currently being developed to describe the magnetic storm 
effects in the inner zone with greater accuracy than presently given 
by AE-5. 
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3. PROGRAM MODEL 


Program MODEL is a Fortran program that enables the user to access 
any of the current trapped radiation models available through the 
National Space Science Data Center (NSSDC) . These models include the 
Inner Zone Electron Model AE-5 for epoch October 1967 described in 
brief in Section 2 of this document and described in detail by Teague 
and Vette (1972), the Outer Zone Electron Model AE-4 for epoch 1967 
given by Singley and Vette (1972) , and a smoothed version of the pro- 
ton models API, AP5 , AP6, and AP7 (originally presented by Vette et al. 
1966-1970) described by Kluge and Lenhart (1971). A matrix storage 
technique originally developed at ESRO (Kluge and Lenhart, 1971) is 
adopted for the containment of these models. A new interpolation scheme 
has been developed at NSSDC and is described in the following sections. 

The matrix storage scheme and interpolation routines are completely 
general, and, as new models become available, these can be easily in- 
corporated into Program MODEL. Work is currently in progress on de- 
veloping new proton and inner zone high-energy electron models. 

A. Program Logic and Restrictions 

Flux versus B/Bo curves are stored in Program MODEL at discrete 
energies and L values using the scheme indicated in Figure 21. Using 
the decadic logarithm of the omnidirectional integral flux, equal 
increments in the ordinate are chosen and the B/B 0 intervals 6(B/Bo)i 
are determined. Each flux versus B/B 0 curve is represented in the 
stored matrix by the variable F 0 , equal to the logarithm of the flux 
at the equator, and the B/Bq intervals <5(B/Bo)i. Using equal incre- 
ments in the ordinate as opposed to the abscissa has the advantage that 
a fixed accuracy is maintained for the flux versus B/B 0 curve even in 
the region of the atmospheric cutoff, where the slope of the curve 
becomes very large. Linear interpolation on the logarithm is used 
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between the grid points defined by <S(B/B 0 )i, and the accuracy of the 
interpolated flux is essentially determined by the ordinate increment. 
Four points per decadic cycle are stored in the present matrix for the 
electron models and two per cycle for the proton model. Flux versus B/B 0 
curves are stored in this manner at a variety of energies and L values. 
Linear interpolation on the logarithm is performed to obtain fluxes at 
intermediate energies and L values. Sufficient energies and L values 
are stored such that an exponential assumption between grid points pro- 
vides sufficient accuracy. This is determined by the radial profiles 
and spectra of the models and is therefore model dependent. The energy 
and L value grid points used for the three models are shown in Figure 
22. Linear logarithmic interpolation between these grid points intro- 
duces less than 10% error in the flux, i.e., considerably less error 
than is presently associated with the models. 

Program MODEL performs the interpolation between grid points in 
the following order: (1) B/B 0 interpolation, (2) L interpolation, and 
(3) energy interpolation. The interpolation scheme adopted for B/B 0 
and L is presented in Figure 23. In this figure, the flux is required 
at some nongrid point P(B/Bo,L) for which the nearest surrounding grid 
L values are Lx and L 2 . A number of rays are drawn from the origin 0, 
taken as B/Bo = 1, log 10 (flux) = 0, to the four grid points surrounding 
P> Ax, A 2 for L = Lx and Bx, B 2 for L = L 2 . The intermediate points 
and A' are determined by linear interpolation on the log xo (flux) 
and B/Bo between points Bx and B 2 and Ax and A 2 , respectively. Further 
linear interpolation is performed to obtain Cx and C 2 at the required 
L value. A final interpolation between Cx and C 2 is performed to 
obtain the correct B/Bo value at point P. In the event that the grid 
L flux-B/Bo distributions cross (as occurs, for instance, in AE-5 at 
low L values and intermediate energies), a number of additional rays 
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are drawn. For nongrid energies, the interpolation scheme described 
above and shown in Figure 23 is used at the two surrounding grid ener- 
gies, and linear interpolation on logio(flux) and E is performed to 
obtain the flux at the correct energy. 

The Kluge and Lenhart (1971) scheme stores the flux using B as 
grid points and interpolates between grid L and E values at constant 
B. The present storage and interpolation scheme has a number of ad- 
vantages over the Kluge and Lenhart method. For L 2 > Li in Figure 23 
at the equator Bo(L 2 ) is less than Bo (Lx) and, at atmospheric cutoff, 
B c (L 2 ) is greater than B c (Li). Thus there are two regions for B < Bo(Li) 
and B > B c (Li) for which flux values can be determined for L = L 2 only. 
Interpolation at constant B/B 0 removes the region B < Bo(Li), but the 
problem remains at the cutoff. With the technique described in the 
previous paragraph, however, the interpolation is performed in a com- 
pletely general fashion without restriction on B/Bo or B. In addition, 
the Kluge and Lenhart scheme is inaccurate for low L values where the 
equatorial B value is quite different from one L grid point to the next. 

In the following paragraphs, a brief description of the main pro- 
gram and subroutines of Program MODEL and their restrictions is given 
in the order in which they are called. 

MAIN 

MAIN performs the I/O function of Program MODEL and offers a 
variety of options to the user for inputting B, L, and E. These options 
are described fully in Section 3B. The variable retrieved by the inter- 
polation subroutines is omnidirectional integral flux with units par- 
ticles/cm 2 -sec. 


11 



The radial profile and spectra for each model have been smoothed at 
grid points, and each model provides fluxes down to 1 particle/cm 2 -sec. 
Smaller fluxes are defined as zero. MAIN is able to determine average 
differential flux in particles/cm 2 -sec-MeV for limited energy ranges. 

An infinite number of grid energies would be required to determine 
smooth point differential fluxes, and practical limits on the energy 
bandwidth result from the finite number of grid energies stored (Figure 
22). These practical limits are determined by imposing the restrictions 
that the resulting differential spectra and radial profiles must remain 
smooth. They are determined to be: 

Protons: E < 1 MeV AE > 250 keV 

1 < E/ (MeV) <20 AE > 1 MeV 

20 < E/(MeV) <50 AE > 5 MeV 

E > 50 MeV AE > 10 MeV 


Inner Zone Electrons: E < 100 keV AE > 50 keV 

100 < E/(keV) < 250 AE > 100 keV 

E > 250 keV AE > 200 keV 

Outer Zone Electrons : same as inner zone electrons except 

E > 4 MeV AE > 100 keV 

An additional restriction is given for AE-5 as L > 1.2 . MAIN 

tests that these conditions are satisfied and disallows smaller energy 
intervals than those shown above. 

Program MAIN supplies the interpolation routines with the particle 
type, a single B and L value, and an array of energies. Multiple B and 
L values are obtained by looping within MAIN. 
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Subroutine TYPE 


Subroutine TYPE is primarily a buffer routine between MAIN and the 
interpolation subroutines that facilitates the incorporation of program 
MODEL into existing programs (Section 3B) . In addition, TYPE determines 
the model to be accessed and converts from logarithm of the flux to 
flux. 

Subroutine TRARA1 

Subroutine TRARA1 determines the grid energies to be retrieved and 
performs the energy interpolation. 

Subroutine TRARA2 

Subroutine TRARA2 determines the grid L values to be retrieved and 
performs the B/Bo-L interpolation shown in Figure 23. 

BLOCK DATA Statements 

The BLOCK DATA statements contain the grid B/Bo, L, and E points 
stored for each model and shown in Figure 22 . Three BLOCK DATA state- 
ments are included, one for each model (see Section 3B for removal of 
unneeded models) . The format of the BLOCK DATA statements is shown 
in Figure 21 for an arbitrary grid energy E. The format is repeated 
for each grid energy. The variables E, L, Fo, and <$(B/B 0 ) are scaled 
such that they can be stored in the BLOCK DATA statements with an 16 
format. The first number N of each grid energy Ej is the total number 
of points in the BLOCK DATA statement corresponding to that energy. 

A general flux versus B/B 0 curve, Jj versus (B/Bo)j, corresponding to 
L - Lj is represented by the number of elements Nj at Lj , the L value 
L j , the decadic logarithm of the equatorial flux, F 0 , and nj B/B 0 in- 
crements where Nj = nj +3. This format is repeated for each grid L 
value. The first two and the last L values stored are end points 
having Fq =0. 
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Subroutine DIFF 


Subroutine DIFF accumulates the average differential flux for 
writing out by MAIN. 

B. Use of Program MODEL 

Versions of the Fortran program MODEL that are suitable for opera- 
tion on IBM 360 series or UNIVAC 1108 computers can be supplied to a 
user. Source deck setups for operation on these machines are shown 
in Figure 24. For operation with source decks, approximate CPU times 
are 6 minutes (IBM 360/75) and (UNIVAC 1108) for 2.5 x 10 1 * points in 
B-L-E space. Reductions in CPU compile time are obtained if Program 
MODEL is executed with object BLOCK DATA statements. Because of varia- 
tions in compiler speed, the actual savings are machine dependent. For 
the Fortran G compiler on the IBM 360/75 at the Goddard Space Flight 
Center, a reduction of a factor of 7 in CPU compile time is obtained 
for 2.5 x lO 4 B-L-E points. In general, however, a factor of 2 re- 
duction may be more typical. Combined object and source deck setups 
are shown in Figure 24. 

Program MODEL offers the user a number of options determined by 
the data cards described in the following paragraphs. 


Card 

Number 

Variable 

Name 

Columns 

Format 

Function 

a 


1-18 

613 



NE 

1-3 

13 

Number of energies . Maximum 
NE = 9 for line printer outpul 
(IPUN = 0) and NE = 5 for cart 
output (IPUN = 1) . Program 
terminates for NE = 0. 


NL 

4-6 

13 

Number of L values. Maximum 
NL = 100 limited by DIMENSION 
statements only. 
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Card Variable 

Number Name Columns Format Function 


MTYPE 7-9 13 Particle type. MTYPE = 2 for 

electrons . MTYPE = 1 for pro- 
tons . 


IDIFF 9-12 13 Determines type of tabular out- 

put. IDIFF = 0 for integral 
flux output, IDIFF = 1 for aver- 
age differential flux, and 
IDIFF = 2 for both. For 
IDIFF = 1 or 2, NDELB (card e) 
is restricted to SO by DIMENSION 
statements . 

IDEF 13-15 13 Determines type of B/Bo range 

used. For IDEF =0, program 
defaults to 25 to 30 linear 
B/Bo increments over the range 
Bo to atmospheric cutoff. 

IDEF = 1 for user input (card e) . 

IPUN 15-18 13 Determines type of output. 

I PUN = 0 for line printer; 

IPUN = 1 for card output (see 
also variable NE) . For 
IPUN =1, output variables are 
L, B, B/Bo, and integral flux 
for each input energy with for- 
mat (F6.2, F8.4, F8.3, 

5(1PE10.3). 

b 1-63 9F7.3 


E F7.3 Energy (MeV) array of length NE 

(card a) . Energies can be in- 
put in any order (see final 
paragraph. Section 3B) . If 
average differential flux is 
required (IDIFF * 1 or 2, 
card a) , this is determined 
in the interval E(I1) to 
E(I + 1) after the E array has 
been sorted into ascending order. 
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Card 

Number 

Variable 

Name 

Columns 

Format 

Function 

c to d 


1-77 

11F7.3 

Number of cards determined by 
NL (card a) . 


XL 

1-7 

8-14 

etc. 

F7.3 

L value array of length NL 
(card a) . L values can be 
input in any order. 

e to f 


1-23 

2E10.3, 13 

These cards are omitted for 
IDEF = 0 (card a). For IDEF = 
one card is required for each L 
value . 


B01 

1-10 

E10.3 

Lower limit of B/B 0 required. 


B0 2 

11-20 

E10.3 

Upper limit of B/Bo required 
for NDELB f 1. For NDELB = 1, 
B02 is the required B/Bo in- 
crement, and the upper limit 
of B/Bo corresponds approxi- 
mately to atmospheric cutoff. 


NDELB 

21-23 

13 

Number of B/B 0 intervals re- 
quired between B01 and B02 for 
NDELB f 1. For NDELB = 0, 
program defaults to NDELB = 20. 

g 


1-18 

613 

As card a. Program terminates 
for NE = 0. 


A summary of the data setup is shown in Figure 25. 


Sample integral and differential flux tables are shown in Figure 
26 obtained with the variables NE = 9, NL = 1, MTYPE = 2, IDIFF = 2, 
IDEF = 0, E = 0.05, 0.1, 0.25, 0.5, 1.00, 1.25, 1.5, 1.75, 2.0, and 
XL = 1.7. In addition to the output shown, a number of messages may 
be printed in association with the model restrictions described in 
Section 3A. If at a given L value no flux is found greater than 1 
particle/cm 2 -sec or cm 2 -sec-MeV, a message is written to that effect. 
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In addition, messages are given if the average differential flux can- 
not be accurately determined because the L value is less than 1.2 or 
because the energy interval input is too small. 

If all BLOCK DATA statements are not required for regular usage 
of Program MODEL, the unneeded models can be removed by removal of (a) 
the appropriate BLOCK DATA statement in which the model is identified 
by the variable NAME, (b) the associated COMMON blocks from MAIN, and 
tc) the associated calls to subroutine TRARA1 made from subroutine 
TYPE and identified by comment cards . The operation of Program MODEL 
remains as described above. 

Program MODEL is designed to be easily incorporated into existing 
programs. A single call to subroutine TYPE is required in the existing 
program : 


CALL TYPE (MTYPE , B, FL, NE, E, FLUX) 

where E, the energy array in MeV, and FLUX, the integral omnidirectional 
flux array returned by TYPE for these energies, must be dimensioned to 
NE in the existing program. In the calling argument for TYPE, MTYPE is 
the particle type as described in card a, and B and FL are, respectively, 
the required magnetic field strength in gauss and the L value in earth 
radii. The variable FLUX has units of particles/cm 2 -sec and has been 
equated to zero for fluxes of less than 1 particle/cm 2 -sec. A single 
additional restriction is imposed upon the user as a consequence of the 
interpolation algorithm used in subroutine TRARA. The energy E must be 
supplied to subroutine TYPE as an ascending array. 
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4. PROGRAM ORP 


The Orbital Radiation Program (ORP) is a Fortran program designed 
to calculate the average geomagnetically trapped radiation accumulated 
by an orbiting vehicle. ORP is a substitute for Program TRECO, pre- 
viously issued by NSSDC (Lucero, 1968), and differs from that program 
in three respects. First, ORP requires B and L coordinates for the sat- 
ellite orbit. Programs for the calculation of the B and L coordinates 
from latitude, longitude, and altitude can be supplied by NSSDC (King, 
1971). Secondly, ORP uses Program MODEL, described in Section 3, for 
determining the particle omnidirectional integral flux along the orbit. 
As noted in Section 3, Program MODEL is general and will be able to 
contain new particle models, as they become available, with only minor 
modification. Finally, ORP does include an orbit generation facility. 

ORP is able to generate the following tabular output. 

Table 1 Intermediate Printout - a point by point printout of the 
omnidirectional integral flux at each point of the orbit 
for a given threshold energy. 

Table 2 L-Band Summary - a summary of the omnidirectional particle 
flux (particles/cm 2 -day) accumulated in arbitrary energy 
and L bands. 

Table 3 Integrated Flux - a summary of the integrated flux accumu- 
lated in arbitrary energy bands. 

Table 4 Intensity Summary - a summary of the omnidirectional parti- 
cle flux accumulated in arbitrary energy and intensity bands . 

Table 5 Peak Flux per Orbit - a table of peak omnidirectional in- 
tegral flux encountered for each revolution for a given 
energy threshold. 

Table 6 Standard Circular Orbits - a summary of omnidirectional 
fluxes (particles /cm 2 -day) to be used only for standard 
circular orbits for four inclinations at a given altitude. 

This information may also be written on tape. 


preceding page 
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Any combination of the above tables may be obtained for a given 
program run with the restriction that the arbitrary energy bands are 
fixed for a given run. Examples of these tables are shown in Figure 
27. 


A. Program Logic 

MAIN 

ORP uses inputted logical controls to determine the types of tab- 
ular output to be presented. A search on the input tape (Section 4B) 
is initiated to locate the first of the orbits needed. ORP loops to 
determine the flux for each point along this orbit. The Intermediate 
Printout table is written out in this loop. At the end of each orbit, 
the summary tables described in Section 4 are written out and the pro- 
gram proceeds to the next orbit or terminates . Each new orbit data 
set must follow the previous set on the tape. For the special case 
of standard circular orbits at 0°, 30°, 60°, and 90° inclination, 

MAIN writes the Standard Circular Orbits table at the end of each 
fourth orbit. This last output table is primarily used for presenting 
exposures along standard orbits for inclusion in model documentation. 
Examples of this table are presented by Singley and Vette (1972) in 
the documentation for the outer zone electron environment AE-4. 

Subroutines TYPE, TRARA1, TRARA2 , and BLOCK DATA Statements 

Hiese subroutines are identical to the ones previously described 
in Section 3A and are the interpolation subroutines and model matrices 
providing omnidirectional integral flux at the B-h-E points supplied 
to TYPE by MAIN. 

Subroutine STORE 

Subroutine STORE accumulates integral fluxes in L bins for the 
L-Band Summary table. The L-bins are specified by the IF statements 
in subroutine STORE. 
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Subroutine FLITAB 


Subroutine FLITAB accumulates integral fluxes in intensity ranges 
for the Intensity Summary table. The intensity ranges are specified 
by DATA statement FLXBIN in subroutine FLITAB. 

Subroutine DECACC 

Subroutine DECACC determines the peak integral omnidirectional 
flux in each revolution for a given threshold energy. This energy is 
user input and is the same energy as used for the Intermediate Printout 
table. For nonequatorial orbits, the south-north crossing of the equa- 
torial plane is used to denote the start of each orbit. For equatorial 
orbits, a local time of zero is used. Subroutine DECACC is not accessed 
if altitude is zero. 

Subroutine DEPRNT 

Subroutine DEPRNT computes the flux accumulated per day in the 
user input energy bands for the L-Band Summary, Integrated Flux, and 
Intensity Summary tables. DEPRNT writes these tables and the Peak 
Flux per Orbit table. 

Subroutine STAND 

Subroutine STAND accumulates the omnidirectional flux per day in 
energy bands for the Standard Circular Orbits table and writes this 
table on the line printer and on tape. 

B. Use of Program ORP 

When Program ORP is run on an IBM 360/75 computer under MVT with 
the Fortran IV H, opt=2 compiler, the run time for four standard (0°, 
30®, 60°, and 90° inclination) orbits of 1440 data points each is ap- 
proximately 2 minutes CPU and 0.5 minutes I/O time using object BLOCK 
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DATA statements. As noted in Section 3B, significant savings in CPU 
time are obtained by using object rather than source decks on a 360/75. 
The execution step used approximately 150 K bytes of storage with the 
two electron models included. A sample deck setup is shown in Figure 
28 for combined object and source decks. 


The orbit information is input to Program ORP using tape input with 
data set reference unit number 10 (see Figure 28) . Each orbit is pre- 
ceded by an alphanumeric header record of up to 76 characters in length 
that describes the orbit and is written out in Tables 1 through 5 
(Section 4A) . The header record is followed by data records giving 
standard geocentric and B and L coordinates for each time. Each record 
contains the following data in E format: longitude, latitude, altitude, 

B, L, and time in hours since start of orbit. An altitude of -100 de- 
notes the end of an orbit. The B and L coordinates are mandatory, but 
the latitude, longitude, and altitude information is optional unless 
the Peak Flux per Orbit table is required. These variables should be 
set to zero or left blank if actual values are not to be supplied. 

The input tape format is as follows: 


Header Record (one per orbit) 
Variable Name Format 


Function 


HEAD 19A4 

Data Record (one per point on orbit) 


Alphanumeric Orbit Description 


Variable Name 

Format 

0RBVAL 

(1) 

E18.8 

ORBVAL 

(2) 

E18.8 

ORBVAL 

(3) 

E18.8 

ORBVAL 

(4) 

E18.8 

ORBVAL 

(5) 

E18.8 

ORBVAL 

(6) 

E18.8 


Function 

Longitude (degrees) 
Latitude (degrees) 
Latitude (km) 

B (gauss) 

L (earth radii) 

Time (hours from start of 
orbit) 
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Trailer Record 


Variable Name Format Function 

ORBVAL (3) E18.8 Set equal to -100 to signal 

end of orbit. Other ORBVAL 
variables are not important 
for this final record. 

The various options of Program ORP may be obtained by use of the data 
deck setup described in the following paragraphs. 

Card Variable 


Number 

Name 

Columns Format 

Function 

a 


1-70 10L1,10A4, 




I3,F6.2, 




10X,I1 



TABCON (1) 

1 

LI 

TABCON(l) = T for Intermediate 
Printout table, F for no table. 

TAB CON (2) 

2 

LI 

TABCON (2) = T for L-Band Summary 
table, F for no table. 

TABCON (3) 

3 

LI 

TABCON (3) = T for Integrated 
Flux table, F for no table. 

TABCON (4) 

4 

LI 

TABCON (4) = T for Intensity 
Summary table, F for no table. 

TABCON (5) 

5 

LI 

TABCON (5) = T for Peak Flux 
per Orbit table, F for no 
table. 

TABCON (6) 

6 

LI 

TABCON (6) = T for Standard 
Circular Orbit table, F for 
no table. 

TABCON (7) 

7 

LI 

TABCON (7) = T for tape output 
of Standard Circular Orbit 
table, F for no table. 

TITLE 

11-50 

10A4 

Alphanumeric array for writ- 


ing at top of first page. 
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Card Variable 


Number 

Name 

Columns 

Format 

Function 


NE 

51-33 

13 

The number of integral 
energy values input. 
Maximum NE = 30. 


ET 

54-59 

F6.2 

The threshold energy (MeV) 
used for the Intermediate 
Printout and Peak Flux per 
Orbit tables. ET may be 
omitted if these tables 
are not required. 


MODEL 

70 

11 

Particle Type. MODEL = 2 
for electrons, 1 for pro- 
tons . Appropriate BLOCK 
DATA statements must be 
included (Section 3). 

b to c 


1-80 

10(F6.2,2X) 

Number of cards determined 
by NE. Ten values per card 

e 


1-6 

F6.2 

Energy threshold array of 
length NE. 



9-14, 

etc. 


Must be in ascending order. 
For tabular output Tables 
2, 3, 4, and 6 (Section 4A) 
the energy intervals E(I+1) 
to E(I) are subject to the 
restrictions given in 
Section 3A under MAIN. 

This function is not per- 
formed automatically by 
ORP. 

d 


1-4 

212 



rOR8(l) 

1-2 


Index number for the first 
orbit required on the input 
tape. The first orbit is 
I0RB(1) = 1. 
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Card 

Number 

d 

Ccont’d) 


A summary 


Variable 

Name Columns 

Format 

Function 


IORB (2) 3-4 

12 

Index number for last 
required on the input 

orbit 

tape. 


If this tape contains a 
single orbit, I0RB(2) may 
be left blank. 

of the data setup is given in Figure 29. 
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Table 1. Omnidirectional Flux Confidence Limits 


Code 

B Range 

L Range 

E_ Range 

Section 

Comment 

1 

> B o 

1 .2-1 *4 

>3 MeV 

4 

Extrapolation on both B dependence and 
spectrum, no data 

2 

" B o 

1.2-1. 4 

>3 MeV 

4 

Extrapolation on spectrum, no data 

4 

> B o 

1.2-1. 7 

>250 keV 

GB 

Possible presence of Starfish electrons 

5 

* B o 

1.9-1. 4 

4-2 MeV 

4 

Magnetic storm effects, single data set, 
B extrapolation 

6 

> B o 

1,7-1. 9 

>500 keV 

4 

Single data set, B extrapolation 

6 

i8 o 

< T.25 

all energies 

5 

L extrapolation 

G 

~ B o 

>1.5 

all energies 

3 

Poor data 

7 

> B o 

1.3 

all energies 

3 

Poor OGO data 

8 

>>B 0 

>2 

all energies 

3 

Poor pitch angle coverage 

10 

* B o 

1.4-1. 9 

< 250 keV 

3 

Agreement between three data sets 
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Table 2. Solar Cycle Parameter Confidence Limits 


Code 

L Range 

e t 

T Range 

Section 

Comment 

3 

<1.8 

250,500 

>22 

6B 

Significant Starfish flux at T=22 re- 
sulting in iteration 

4 

>1.9 

all 

all 

3 

Poor QG0 data at high L values 

4 

all 

40 

all 

6A 

Small Rf values; data standard deviation 
becomes significant. 

5 

<1.4 

250,500 

all 

6B 

Hardening of spectrum; assumed constancy 
of j (>690 ) term in equation 16 

5 

>1.8 

500 

all 

6A 

Storm effects term In equation 17 becomes 
signi ficant 

7 

1. 6-2.0 

250,500 

all 

6A 

Two data sets available {i.e., QGO and 
1963-38C) 


Note that these confidence codes are low because integral flux values of Ry are determined from the 
0G0 data using an approximate expression, and B independence has been assumed. Further, if Ry is 
used to extrapolate beyond the epoch of October 1967, as described in Section 7, the above confi-' 
dence codes will be reduced because of asymmetries in the solar cycle. 




Table 3. Cutoff Times for Starfish Electrons 



p = 0.5 


P = 

0.25 


L 

Ch 2 

Ch 3 

Ch 4 Ch 5 

Ch 2 

Ch 3 

Ch 4 

Ch 5 

1.4 

2/65 (31 ) ±3 

7/66 (48) ±6 

3/68 (68) ±8 1/68 (66) ±7 

8/65 (37) 

11/66 (52) 

7/69 (84) 

12/68 (77) 

1.5 

3/65 (32) ±3 

7/66 (48) ±6 

9/67 (62) ±7 4/67 (57) ±6 

9/65 (38) 

11/66 (52) 

10/68 (75) 

3/68 (68) 

1.6 

3/65 (32) ±2 

7/66 (48) ±4 

9/66 (50) ±2 3/66 (44) ±2 

10/65 (39) 

11/66 (52) 

8/67 (61) 

2/67 (55) 

1.7 

4/65 (33) ±2 

2/66 (43) ±3 

2/66 (43) ±2 -10/65 (39) ±2 

12/65 (41) 

8/66 (49) 

8/66 (49) 

-4/66 (45) 

1.8 

3/65 (32) ±2 

2/66 (43) ±3 

ND <12/64 (<29) 

10/65 (39) 

6/66 (47) 

ND 

<12/64 (< 29) 

1.9 

1/65 130). ±2' 

11/65 (40) ±3 

ND <12/64 (<29) 

9/65 (38) 

4/66 (45) 

ND 

<12/64 (<29) 

2.0 

1/65 (30) ±2 

ND 

ND <12/64 (<29) 

7/65 (36) 

ND 

ND 

<12/64 (< 29) 

2.2 

10/64 (27) ±2 

ND 

ND <12/64 (< 29 ) 

2/65 (31) 

ND 

ND 

<12/64 (< 29) 


p * Vj 


ND denotes no data because of magnetic storm effects or no measurements. 


Figures in parentheses represent months from Starfish injection. 


Channel 

2 

3 

4 

5 


Energy Range (keV) 

133-292 

292-690 

690-1970 

1970-4740 




Figures 1-7 . These computer-generated plots present carpet plots 
of the AE-5 omnidirectional flux as functions of B and L for threshold 
energies E T = 0.04, 0.1, 0.25, 0.5, 1.0, 2.0, and 4.0 MeV. A descrip- 
tion of the use of these carpet plots is given in Appendix A of NSSDC 
72-10 (Teague and Vette, 1972). In general, lines of constant B are 
presented in 0.02-gauss increments from the equator to 0,28 gauss, and 
lines of constant L are presented in increments of 0.05 earth radii for 
1.2 < L < 1.6 and increments of 0.1 earth radii for 1.6 < L < 2.4. In 
some cases, lines are omitted for clarity. For the energies 1.0, 2.0, 
and 4.0 MeV, the plots are subdivided into two L ranges because of 
steep gradients in the radial profiles at these energies. In each 
figure the ordinate scale increments are shown as error bars on the 
left-hand side of the plot and the abscissa scale is shown as powers 
of ten. 
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FIGURE 1 

AES OMNIDIRECTIONAL INTEGRAL FLUX, ET=40 KEV 
EPOCH OCTOBER 1967 
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AES OMNIDIRECTIONAL INTEGRAL FLUX, ET=100 KEV 
EPOCH OCTOBER 1967 
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FIGURE 3 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET=250 KEV 
EPOCH OCTOBER 1967 
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FIGURE 4 

AES OMNIDIRECTIONAL INTEGRAL FLUX, ET=500 KEV 
EPOCH OCTOBER 1967 
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FIGURE 5_ 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET= 1 MEV 
EPOCH OCTOBER 1967 
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FIGURE S CONT 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET= 1 MEV 
EPOCH OCTOBER 1967 
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AE5 OMNIDIRECTIONAL INTEGRAL FLUX* ET=2 MEV 
EPOCH OCTOBER 1967 
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CONT 


AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET= 2 MEV 
EPOCH OCTOBER 1967 
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FIGURE 7 

AES OMNIDIRECTIONAL INTEGRAL FLUX, ET=4 MEV 
EPOCH OCTOBER 1967 
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CONT 


AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET=4 MEV 
EPOCH OCTOBER 1967 
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Figure 13. AE-5 R-X Flux Map (electrons/cm 2 -sec) 
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FIGURE 14 


A E 5 OMNI-DIRECTIONAL INTEGRAL FLUX 
GREATER THAN 40 KEV 
EPOCH OCTOBER 1887 
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Figures 14-16. Three-Dimensional Flux Maps 
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4 3 2 1 

r t 

Figure 17. Integral Flux Solar Cycle Ratios Ry, E > 40 keV 
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Figure 18. Integral Flux Solar Cycle Ratios Rj, E > 100 keV 
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Integral Flux Solar Cycle Ratios Rj, E > 250 keV 
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Protons 


Energies: 0.375, 0.78, 4.1, 8.0, 16.0, 50.0, 100.0 MeV 

L values: 1.2 by 0.1 increments to 6.6 

Inner Zone Electrons 

Energies: 0.04, 0.1, 0.25 by 0.25 increments to 2.0, 

2.0 by 0.5 increments to 4.5 MeV 
L values: 1.2 by 0.05 increments to 1.5, 1.5 by 0.1 

increments to 2.0, 2.0 by 0.2 increments to 2.8 

Outer Zone Electrons 

Energies: 0.04, 0.1, 0.3, 0.5, 1.0, 2.0, 2.5, 3.0, 3.5, 

4.0, 4.1, 4.25, 4.35, 4.5, 4.65, 4.85 MeV 
L values: 2.8 by 0.2 increments to 4.0, 4.0 by 0.5 increments 

to 6.0, 6.6, 7.0 by 1.0 increments to 11.0 


Figure 22. Model Grid Points 



LOG -|o (flux) 



B/B 0 


r lyurc 


23 . 


n /n 

d/d 0 


aHu 


T 

in 


te 


rpCu atl un 


60 


IBM 360 Series 


Univac 1108 


Job Card 

// EXEC F0RTRANG,PARM= ' ID, MAP, 
XREF' ,REGI0N=25QK + 

//SOURCE, SYS IN DD * 

Source Deck 
/* 

//STEPG EXEC LINKG0,REGI0N=160K 

//LINK. OBJECT DD * 

Object Deck if used 

/* if Object Deck used 

//GO.SYSUDUMP DD SYS0UT=A 

//GO.GSFCDUMP DD SYS0UT=A 
CGSFC only) 

//G0.DATA5 DD * 

Data Deck (.Figure 25) 

/* 

// 


Job Card 

0 F0R,SIA .MAIN..MAIN/R 
Main Program 

0 F0R,SIA .SUB1 ,.SUB1/R 
Subroutine Type 

0 F0R,SIA .SUB2..SUB2/R 
Subroutine DIFF 

0 F0R,SIA ,SUB3,.SUB3/R 
Subroutine TRARA1 

@ F0R,SIA .SUB4,.SUB4/R 
Function TRARA2 

0 F0R,SIA .SUB8,.SUB8/R 
AE5 Block Data 

6 F0R,SIA .SUB9 , .SUB9/R 
AE4 Block Data 

0 FOR.SIA .S010..S010/R 

0 MAP, I .MAIN/R, .MJTP/A 

0 XQT .MJTP/A 
Data Deck 

0 FIN 


+Source block data statements only. 


Figure 24. MODEL Deck Setup 
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Figure 25. MODEL Data Deck Setup 
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Figure 26. Sample MODEL Output 
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Figure 26 (continued) 
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Figure 27. Sample ORP Output 
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Figure 27 (continued) 
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Figure 27 (continued) 



PEAK FLUX PER GR0I T TABLE 

CIRCULAR ORBIT* PER ft AP * 300 NX. FIELD N00 - HftC-120. INCL »30 DEG 

NODELS USED * AE4* AES 
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0.0 

217.17 

0.76 

555. 6 

9.SA333 

0.0 
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Figure 27 (continued) 



four orbit Integrated flux summary 


MODEL USEO-AEAt AES 


OREIT ALTITUDE** 300 N MI 


TOTAL TIME.* 34* HOURS 


TIME INTERVAL*. 60. SECONOS 


ENERGY 

MEV 


ORBITAL FLUX 
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Figure 27 (continued) 



Job Card 

// EXEC FORTRANH ,PARM= 1 0PT=2 , I D .MAP ,XREF ' ,REGI0N=450K + or 375K X 
//SOURCE. SYSIN DD * 

Source Deck 
/* 

// EXEC LINKGO.REGION.GO = 160K 
//LINK. OBJECT DD * 

Object Deck 
/* 

//GO . FTTOFOOl DD (input tape information) 

//GO.FT13F001 DD (output tape information) 

//GO .SYSUDUMP DD SYSOUT=A 

//GO .GSFCDUMP DD SYSOUT=A (GSFC only) 

//G0.DATA5 DD * 

Data Deck (Figure 29) 

/* 

// 


+Source block data statements. 
^Object block data statements. 


Figure 28, ORP Deck Setup 
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Figure 29. ORP Data Deck Setup 
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